r/n- 


NASA Technical Memorandum 81856 


NASA-TM-81856 19810019011 


COMPUTER PROGRAM FOR NONLINEAR STATIC STRESS ANALYSIS 
OF SHUTTLE THERMAL PROTECTION SYSTEM - USER'S MANUAL 


Gary L. Giles and 
Maria Valias 



HOI » S '*"' 


Hay 1981 

JUN 2 2 1981 

CENTF.1^ 
LIBt^AaV, WA3/\ 

vippiiJ'i'wX 



t\i/^ 

National Aeronautics and 
Space Administration 

Langley Research Center 

Hampton, Virginia 23665 




CONTENTS 


Page 

SUMMARY 1 

INTRODUCTION 1 

ANALYSIS PROCEDURE • 2 

DESCRIPTION OF COMPUTER PROGRAM 4 

SUMMARY OF COMPUTER PROGRAM INPUT 5 

DESCRIPTION OF INPUT COMMANDS 7 

Case Identification 7 

Command: "TITLE" 7 

Command: "CASE ID" 7 

Command: "COMMENT" 8 

Geometry Description 8 

Command: "NODE" 8 

Command: "CORNER" 9 

Command: "NORMAL" 10 

Command: "POINT" H 

Command: "PLANE-LINE" 12 

Command: "LOCAL" 13 

Command: "BLOCK" 13 

Command: "WEDGE" 14 

Command: "TETRA" 18 

Command: "FACE" 16 

Command: "CG" 17 

Command: "VIEW" 17 

i 






Page 

Tile Loads 18 

Command: "FORCE" 18 

Command: "MOMENT" 18 

Command: "PRESSURE" 18 

Command: "GLOAD" 19 

Nonlinear Analysis Data. , 19 

Command: "MATERIAL" 20 

Command: "IMPERFECT" 24 

Command: "PAD" 29 

Command: "CONTROL" 32 

Command: "NONSTRESS" 34 

EXAMPLE ANALYSIS 35 

CONCLUDING REMARKS 48 

APPENDIX A - MATERIAL CURVE DATA FILE 49 

REFERENCES 52 

FIGURES 53 


i i i 




SUMMARY 


User documentation is presented for a computer program which considers the 
nonlinear properties of the strain isolator pad (SIP) in the static stress 
analysis of the shuttle thermal protection system (TPS). This program is 
generalized to handle an arbitrary SIP footprint including cutouts for instru- 
mentation and filler bar. Multiple SIP surfaces can be defined to model tiles 
in unique locations such as leading edges, intersections, and penetrations. 

The nonlinearity of the SIP is characterized by experimental stress-displacement 
data for both normal and shear behavior. Stresses in the SIP are calculated 
using a Newton iteration procedure to determine the six rigid body displace- 
ments of the tile which develop reaction forces in the SIP to equilibrate the 
externally applied loads. This user documentation gives an overview of the 
analysis capabilities, a detailed description of required input data and an 
example to illustrate use of the program. 

INTRODUCTION 

The metal primary structure on the external surface of the space shuttle 
orbiter is protected from the thermal environment by an array of RSI (reusable 
surface insulation) tiles. These tiles are not attached directly to the shuttle 
surface or substrate but are mounted on an intermediate strain isolator pad as 
shown in figure 1(a). This pad prevents strain in the shuttle structure pro- 
duced by thermal or aerodynamic loading from cracking the brittle tile material. 
In the baseline TPS configuration, the tiles have uniform properties through the 
thickness (undensified tiles). During TPS tests with undensified tiles under 
externally applied loads, failure was found to occur in the tile adjacent to 
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the plane where it was bonded to the SIP. Therefore in this case, the failure 
stress can be taken as the through- the-thickness or normal stress in the SIP. 

As indicated in reference 1, the SIP has nonlinear stress-displacement 
properties. 

The present report contains user documentation for a computer program which 
considers these nonlinear properties in a static analysis for calculating 
stresses in the SIP. The tile attached to the outside of the SIP is considered 
rigid and can have an arbitrary shape built up with block elements as shown in 
figure 1(b). This program is generalized to handle an arbitrary SIP footprint 
including cutouts for instrumentation (fig. 1(c)). Also, filler bar material 
which is usually located under the tile around the perimeter of the SIP but not 
attached to the tile can be considered. Multiple SIP and filler bar areas can 
be defined to model tiles with nonplanar surfaces in unique locations such as 
leading edges, intersections, and penetrations. 

A brief description of the procedure used to calculate the SIP stresses is 
presented in the next section. A detailed description of the required input 
data and an example to illustrate the input and corresponding calculated results 
are given in the remainder of this report. 

ANALYSIS PROCEDURE 

The analysis procedure used in the computer program is a generalization and 
extension of the basic method described in reference 2. The SIP material has 
nonlinear stress-displacement relationships in both the through-the-thickness or 
normal and shear directions. These nonlinear curves are obtained from material 
tests and are represented in the program as tables of stress-displacement 
points. The tile is assumed to be a rigid body with six degrees of freedom 
(three translations and three rotations) at its center of gravity. External 
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loads are applied to the tile as concentrated forces, pressures, and/or stati- 
cally equivalent accelerations (g-load). The substrate under the SIP can be 
given a prescribed shape to represent mismatch from the manufacturing process, 
warpage of the tile, and/or deformations of the external surface of the struc- 
ture such as those caused by buckling. 

A Newton iteration procedure is used to calculate the displacements and 
rotations of the rigid tile for which the reaction forces from stresses in the 
SIP and filler bar material are in static equilibrium with the applied loads. 

The iteration is started using an assumed set of tile displacements and rota- 
tions. This position of the tile relative to the substrate is then used to 
calculate normal and shear displacements of the SIP material. The corresponding 
SIP stresses are obtained from the nonlinear stress-displacement curves. 

Stresses over each SIP surface are numerically integrated and the resulting 
forces and moments are summed at the center of gravity of the tile. The SIP 
reactions are compared to the external loads and the unbalanced or residual 
forces are calculated. During the numerical integration process, a Jacobian 
matrix is formed whose terms are the derivatives or changes in forces and 
moments produced by changes in the tile displacements and rotations. This 
Jacobian matrix and set of unbalanced forces is used to calculate an updated 
set of displacements and rotations by Newton's method. The iteration is con- 
tinued until the SIP forces and external forces are equilibrated. Stresses in 
the SIP are then obtained from the stress-displacement material curves for the 
resulting values of local displacements over each SIP surface. 

To provide capability to analyze all tile configurations the SIP footprints 
are defined with an arbitrary boundary of linear segments including cutouts for 
instrumentation or other penetrations. A general numerical integration method 
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is provided to accommodate SIP areas with complex boundaries. In this method, 
the SIP and/or filler bar surfaces are divided into triangular regions. Within 
each region, the integration process assumes a linear distribution of the inte- 
grand over subtriangles which are defined using area coordinates as discussed in 
reference 3. 


DESCRIPTION OF COMPUTER PROGRAM 

The nonlinear SIP stress analysis program is operational on both the CDC 
NOS 1.3 and the CDC NOS/BE operating systems. The program is organized as a 
primary overlay (0,0) and two secondary overlays (1,0) and (2,0). The primary 
(0,0) and the (1,0) secondary overlays are used for input and generation of 
required geometry and applied load information. The majority of the subroutines 
in these two overlays were written by A1 Dobrowski at Rockwell International for 
use in a linear SIP stress analysis program, reference 4. These subroutines 
were adapted with slight extensions and modifications where necessary for use 
in the nonlinear SIP analysis program. Overlay (2,0) contains the NASA-LaRC 
developed subroutines for implementing the nonlinear analysis procedure 
described in the previous section. The analysis program uses the SPAR data 
handling utilities, reference 5, to access the experimental stress-displacement 
curves which are maintained on a separate data file. The program generates a 
complete file of calculated data on TAPE6 which is usually output on a high 
speed printer. Additionally, a synopsis of stress distributions over each SIP 
pad is output on TAPEl for viewing on an interactive terminal upon completion 
of the analysis run. A single line containing the maximum stress value from 
each analysis with some identification data is output on TAPE2 to form a table 
for runs with multiple tiles and/or multiple load conditions. 
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SUMMARY OF COMPUTER PROGRAM INPUT 

The input data to the program is generally in list directed input form on 
80-column card images as described in reference 6. In list directed form, the 
input data are free-field (can be located anywhere on the card image) and values 
are separated by one or more blanks, or by a comma or a slash, either of which 
may be preceded or followed by any number of blanks. The input deck contains 
commands beginning with a special designator word to indicate which of the 
various types of input data will follow. These commands are described herein 
in the general order that they would occur in an input deck. 

Three commands are available for providing descriptive data to the program. 
These commands are: 

"TITLE" 

"CASEID" 

"COMMENT" 

Geometry input data are designated by the next group of commands. 

Commands used to define node points are: 

"NODE" 

"CORNER" 

"NORMAL" 

"POINT" 

"PLANE-LINE" 

"LOCAL" 

Commands used for calculation of weight and center-of-gravity location of 
the tile are: 

"BLOCK" 

"WEDGE" 
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"TETRA" 


"FACE" 

"CG" 

A command to select a printer generated plot of geometry data is: 

"VIEW" 

The externally applied tile loads are specified by the commands: 

"FORCE" 

"MOMENT" 

"PRESSURE" 

"GLOAD" 

The nonlinear stress analysis requires data given by the following commands: 
"MATERIAL" 

"IMPERFECT" 

"PAD" 

"CONTROL" 

"NONSTRESS" 

Generally, an input card image contains a command designator word followed 
by a list of data values for that command. The designator word must be enclosed 
by quotation marks ("). The list of data values can be continued on more than 
one card image. A slash (/) is needed to terminate a list of data values whose 
length is case dependent such as the number of nodes used to define a SIP 
boundary. The data must be input in a logical sequence or order. For example, 
before a node point can be referenced in a command, it must have been previously 
defined. Hence, the "NODE" data must precede the "CORNER" data. 
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DESCRIPTION OF INPUT COMMANDS 

A detailed description of the input commands is given in this section in 
the general order that they would logically occur in an input deck. 

Case Identification 

Command: "TITLE" 

Purpose: To input a card containing an identification title which is 

printed as a heading on various pages of an analysis output 
listing. 

Format: "TITLE" 

THIS IS AN EXAMPLE TITLE CARD 

Description: The card following the "TITLE" command contains alphanumeric 

information in the first 72 columns. 

Command: "CASEID" 

Purpose: To input a card containing case identification information to 

be output on a file of calculated stress values. These 
values are required by the program for identification pur- 
poses on a special output file, TAPE2, but have no other 
effect on the program. 

Format: "CASEID" NPART XfiACH Q ALPHA 

Description: This is a single card containing the following values: 

NPART - Integer tile part number 
XMACH - Mach number 
Q - Dynamic pressure 
ALPHA - Angle of attack 
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Command: "COMMENT" 


Purpose: 


Format: 

Description: 


Command: "NODE" 

Purpose: 

Format: 

Description: 


To allow users to insert comments directly into the input 
deck. These comments are printed in the output listing as 
they are encountered and have no other effect on the opera- 
tion of the program. 

"COMMENT" N 

N comment cards 

The integer N specifies the number of comment cards, con- 
taining alphanumeric information, which are to be input 
following the "COMMENT" command. 

Geometry Description 

To input nodal coordinates used for the geometric definition 
of the TPS model . 

"NODE" N X Y Z 

This command contains the following values: 

N - Integer node number assigned by user. Maximum value 
which can be used is 200. 

X,Y,Z - The x,y,z coordinates of a node point, 
respectively. 
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Command: "CORNER" 


Purpose: 


Format: 


Description: 


To establish a node point in the corner of a plane through 
previously defined node points. This command is often 
used to define SIP boundary nodes on the inner surface of 
the tile. 

"CORNER" I J K L XI X2 



This command is used to calculate the coordinates of a 
node, L, at given perpendicular distances XI and X2 from 
the edges of a corner I-J-K. 
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Command: "NORMAL" 

Purpose: 

Format: 


Description: 


To establish a node point at a given distance normal to a 
plane. 

"NORMAL" I J K L X 


L 



The coordinates of a point, L, are calculated at a perpen- 
dicular distance, X, from plane I-J-K. The direction of 
J-L is given by the right-hand rule: J-I x J-K. 
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Command: "POINT" 


Purpose: 


Format: 


Description: 


To establish the position of a node point relative to 
previously defined nodes. 

"POINT" I J K L X 



The coordinates of a node, L, are calculated at a distance, 
X, from node, K, in the direction I-J. 
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Command: "PLANE-LINE" 


Purpose: 


Format: 


Description: 


To establish a node point at the intersection of a 
previously defined plane and line. 

"PLANE-LINE" I J K L M N 


K 



I 


The coordinates of a node, N, are calculated at the inter- 


section of plane I-J-K and line L-M. 



Command: "LOCAL" 


Purpose: 

Format: 

Description: 


Command: "BLOCK" 

Purpose: 

Format: 


Description: 


To print the coordinates of all node points with respect to 
a user specified local coordinate system. 

"LOCAL" L M N 

This command contains the following values: 

L - A node defining the origin of the local coordinate system. 

M - A node defining the direction of the local y-axis, L-M. 

N - A node located in the local x-y plane. 

To input the node numbers at the corners, and the density, of 
a block element used for modeling a tile. The sides of the 
block need not be parallel but they must be planar. 

"BLOCK" IJKLMNOP/ DENSITY 


N 0 



The eight integer values I through P are previously defined 
node numbers at the corners of a block listed in the order 
as shown in the sketch. DENSITY is the weight per unit 
volume of the block. 
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Command: "WEDGE" 


Purpose: 


Format: 


To input the node numbers at the corners, and the density, 
of a wedge element used for modeling a tile. Multiple 
BLOCK, TETRA, and WEDGE commands can be used to build up 
irregular shaped tiles. 

"WEDGE" I J K L M N / DENSITY 



Description: The six integer values I through N are previously defined 

node numbers at the corners of the wedge. The triangular 
faces need not be parallel but the quadrilateral faces 
must be planar. DENSITY is given as weight per unit 
volume. 
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Command: "TETRA 


Purpose: 


Format: 


Description: 


To input the node numbers at the corners, and the density, 
of a tetrahedron element used for modeling a tile. 
"TETRA" I J K L / DENSITY 


L 



The four integer values I through L are previously defined 
node numbers at the corners of the tetrahedron. 

DENSITY is given as weight per unit volume. 
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Command: "FACE" 


Purpose: 


Format: 


Descri ption : 


To input the node numbers at the corners, and the density, 
of a face element used for modeling the coating on a tile. 
Each face element must be planar. 

"FACE" I J K . . . / THICKNESS DENSITY CNODE 



Up to 50 previously defined node numbers may be used to 
specify the corners of a face element. 

THICKNESS gives the thickness of the face. 

DENSITY is given as weight per unit volume. 

CNODE is a node number whose coordinates are established 
at the centroid of the face by the face command. 
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Command: "CG" 


Purpose: 

Format: 

Description: 


Command: "VIEW" 

Purpose: 

Format: 

Description: 


To establish a node at the center of gravity of the tile and 
cause this c.g. location and tile weight to be printed. 

"CG" NODE 

This command is used after the tile has been defined as a 
composite of BLOCK, TETRA, WEDGE, and FACE commands. 

The node, NODE, is established at the center of gravity of 
the tile. 

To generate a plot of the tile model for use in checking 
input geometry. This plot is produced by the printer in 
the output listing. 

"VIEW" L M N 

The viewing plane for the plot is established by the 
previously defined nodes L, M, and N. 

L - A node defining the origin of the viewing, x-y, plane. 

M - A node defining the direction of the y-axis, L-M. 

N - A node located in the x-y plane. 
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Tile Loads 


Loading is applied to a tile by combinations of FORCE, MOMENT, PRESSURE, 
and GLOAD commands. A load case number is specified on each of these commands 
and all load contributions for commands with the same load case number are 
summed to form a total load case. Up to fifty load cases can be defined. 
Command: "FORCE" 

Purpose: To input a concentrated force at an existing node point. 

Format: "FORCE" LDCASE NODE FX FY FZ 

Description: A concentrated force with components FX, FY, and FZ along 

the X,Y,Z axes of the global coordinate system is applied 
at node point, NODE, and used in load case, LDCASE. 

Command: "MOMENT" 

Purpose: To input a concentrated moment at an existing node point. 

Format: "MOMENT" LDCASE NODE MX MY MZ 

Description: A concentrated moment with components MX, MY, and MZ about 

the X,Y,Z axes of the global coordinate system is applied 
at node point, NODE, and is used in load case, LDCASE. 

Command: "PRESSURE" 

Purpose: To input a pressure loading on a tile surface. 

Format: "PRESSURE" LDCASE L PL M PM N PN 

Description: A linear pressure distribution is defined over a triangular 

surface with corner node points L, M, and N. The pressure 
over the triangle varies linearly between pressures PL, 

PM, and PN which are specified at the corners with the 
direction of positive pressure given by L-M x L-N. This 
pressure loading is used in load case, LDCASE. 
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Command "GLOAD 


Purpose: To input an acceleration or g-load factors at the tile 

center of gravity. 

Format: "GLOAD" LDCASE GX GY GZ 

Description: The g-loading is specified as components GX, GY, and GZ in 

the global coordinate system. 

Nonlinear Analysis Data 

The nonlinear analysis procedure requires the geometrical definition of the 
SIPS and filler bars which are defined by "PAD" commands and their nonlinear 
material curves defined by "MATERIAL" commands. Substrate deformation shapes 
are defined using "IMPERFECT" commands and the load case combination to be used 
for a particular analysis is given on the "NONSTRESS" command which causes the 
nonlinear analysis to be performed. A "CONTROL" command can precede the 
"NONSTRESS" command to reset parameters which govern execution of the analysis 
program. An index value must be assigned by the user to each material and each 
imperfection. These index values are referenced on the "PAD" commands to 
define the material and imperfections to be associated with each pad. 
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Command: "MATERIAL" 


Purpose: To input the stiffness characteristics of the SIP and/or 

filler bar material. This input is in the form of curves 
which relate stresses to displacements or strains in the 
material. Curves for the SIP material must be specified 
for both normal and shear behavior. Filler bars cannot 
transmit any tensile or shear forces since they are not 
attached to the tile. Therefore, forces are only produced 
by the compression portion of the normal curve with the 
tensile portion represented by a horizontal line at zero 
stress and a shear curve is not required. All material 
curves are extrapolated as straight lines through the last 
two points on each extremity of the curve for cases with 
displacements beyond the bounds of the specified curve. 

Up to twenty different material curves can be defined. 
Three format options are available for input of material 
characteristics: 

(1) Input of a table of points on a stress-displacement 

curve 

(2) Select a material curve from a previously defined data 

file 

(3) Input of linear material characteristics. 

Format (1): "MATERIAL" INDEX "TABLE" T NPTS 

Description: This form of the "MATERIAL" command is used for direct input 

of a table of points to define either a normal or shear 
stress-displacement curve. 
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INDEX - An integer 1 -> 20 used to identify a particular 
material curve. 

"TABLE" - Indicates a table of points follows. 

T - Is the thickness of the material. 

NPTS - The number of points on the curve which follow 

this "MATERIAL" command. One point is given per 
row in the table with the displacement value 
given first and the stress value given second. 
The definition of the curve should start with the point 
having the largest negative displacement value and proceed 
to largest positive displacement. 

Example: 

"MATERIAL" 5 "TABLE" 0.160 6 



-.04 

-15.0 



-.02475 

- 6.0 



-.005 

- 1.0 



0.0 

0.0 



.02178 

5.0 



.04410 

15.0 


Format (2): 

"MATERIAL" INDEX "NORMAL" "Al" "A2" 

N3 



or 




"SHEAR" 


Description: 

This form of the 

"MATERIAL" command is used to 

select a 


stress-displacement curve from a previously defined data 
file. The curves are given a unique four word identifier. 
The procedure used to establish or add curves to the 
required data file is given in Appendix A. 
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INDEX - An integer 1 20 used to identify a par- 

ticular material . This index is referred 
to on the "PAD" command to designate its 
material. 

"NORMAL" 

or - "NORMAL" indicates the material is to be 

"SHEAR" 

used for normal stress calculation. 
"SHEAR" indicates the material is to be 
used for shear stress calculation. 

"Al" and "A2" - Are the first and second words in the 
identifier name containing up to four 
characters and beginning with an alpha- 
numeric character. 

N3 and N4 - Are the third and fourth words in the 

identifier name containing up to four 
integer characters. 


Example curve 

names : 



"FBN" 

"N16C" 

1000 

0 

"SIPN" 

"N16T" 

1580 

100 

"SIPN" 

"R09V" 

0 

0 
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Format (3): 
Description: 


"MATERIAL" INDEX "LINEAR" SLOPE • T 
This form of the "MATERIAL" command is used to define a 
linear stress- strain relation for use in analysis. The 
line defined by this command passes through the origin of 
a stress-strain plot and has the slope which is specified. 
INDEX - An integer 1 20 used to identify a particular 

material curve. 

"LINEAR" - Indicates the material has a linear stress- 
strain relationship. 

SLOPE - Is the slope of the stress-strain curve 
(e.g. , E or G) . 

T - Is the thickness of the SIP material 

(e.g. , .09 or .16) . 
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Command: "IMPERFECr' 


Purpose: 


Format: 

Description: 


To define imperfections such as mismatch, warpage, or 
substrate deflection to be used in analysis. If no 
"IMPERFECT" commands are input, the inner tile surface 
and the substrate surface are taken to be flat. 

"IMPERFECT" INDEX ITYPE I J AIMP GAMY GAMX 

Imperfections are applied to only those SIPS and/or filler 
bars in which they are referenced on the "PAD" commands. 

Up to twenty different imperfections may be defined with 
the imperfection used in a particular analysis being the 
sum of the imperfections specified on the "NONSTRESS" 
command. The command contains the following values: 

INDEX - An integer (1 20) used to give an identification 

number to each imperfection. 

ITYPE - An integer (1 ^6) used to specify an imperfection 
type. The available types are: 

= 1 For cylindrical imperfection 
= 2 For spherical imperfection 
= 3 For single cosine 
= 4 For double cosine 
= 5 For step function 
= 6 For ramp function 

I - Locates the origin of the local coordinate system 
of the imperfection. 

J - Defines the direction of the local x-axis, I-J. 
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AIMP - Gives the imperfection amplitude; a positive value 
indicates a direction into the SIP and a negative 
value is away from the SIP. 

GAMY, GAMX - Wavelength parameters associated with the 
imperfections. 

The meaning of these values that are applicable for 
defining each type of imperfection is given on the 
following figures. Some values are not required in the 
definition of some types of imperfections, but they must 
be input as dummy values. 
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Command: "PAD" 


Purpose: 


Format: 


Description: 


To input information needed to describe a pad of SIP or 
filler bar. Multiple pads can be defined to model tiles 
having a nonplanar attachment surface or to provide for a 
SIP with a variation in material properties. The variation 
could be modeled as a series of discrete pads. 

"PAD" INDEX "SIP" MATN MATS I J K .../IMPl IMP2 .../ 
or 
"FB" 

Up to twenty pads can be defined and each command contains 
the following values: 

INDEX - An integer 1 ^ 20 used to identify a particular 

pad for reference in the output listing. 

"SIP" - "SIP" indicates the pad is SIP material, 
or 

"FB" "FB" indicates the pad is filler bar material. 

MATN - An integer 1 ->■ 20 referring to an entry defined by 
a "MATERIAL" command. The corresponding material 
properties are used for normal stress calcula- 
tions on this pad. 

MATS - Similar to MATN, but refers to material properties 
for shear stress calculations. For filler bars, 
this entry should be zero. 

I J K.../ - A list of nodes which define the boundary of 

the pad. The first three nodes in the 
list set up a local coordinate system for 
the pad. The origin is located at the 
first node and the Y-axis passes through 
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the second node. The third node deter- 
mines the direction of the X-axis (i.e., 
which side of the Y-axis it will be). 

The nodes around the outside of a pad should 
be input in a clockwise direction so that 
the local Z-axis points into the tile. 

For a pad with a hole, such as a filler 
bar, the inside boundary should be input 
in a counterclockwise direction to subtract 
out the contribution of the interior hole. 
Up to two hundred nodes can be input in 
this list which must be terminated with a 
slash. 

An example of numbering a filler bar follows: 



1 4 


Numbering: 

1 234 158765 / 
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IMP! IMP2.../ - A list of integers 1 ^ 20 referring to 

entries defined by "IMPERFECT" commands. 
These imperfections will be the ones 
applied to this particular pad if ref- 
erenced on a "NONSTRESS" command. If no 
imperfections are to be considered for 
the pad, a zero must be input with the 
list terminated with a slash. 
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Command: "CONTROL 


Purpose: The "CONTROL" command is used to define parameters which 

control the nonlinear stress analysis procedure and the 
resulting printed output. 

Format: "CONTROL" IPRT NDIV NSTEPS MESH TEST NCON 


Default 

values 


5 10 .01 0 


If a "CONTROL" command is not input, the default values shown 


above will be used. 


Description: An iterative solution procedure is used to determine the tile 

displacements which give equilibrium between the applied 
loads and reaction of the integrated material stresses in 


the pads. The accuracy of solution and required computa- 
tional time are controlled by the values of NSTEPS, MESH, 
and TEST. The quantity of printed output is controlled by 
the values of IPRT, NDIV, and NCON. 

During the solution procedure, the total applied tile load, 
which is specified as the combination of load cases on a 
"NONSTRESS" command, can be applied as a sequence of equal 
load increments or steps. A converged solution is obtained 
at each load step. The following values control this solu- 
tion process: 

IPRT - = 0 Print iteration history and resulting stresses 
at the final load step which is the total 
applied load. 

= 1 Print iteration history at each of NSTEPS load 
steps and print stresses for only the final 
load step. 



= 2 Print both iteration history and resulting 
stresses at each load step. 

NDIV - Specifies the number of points along edges of pads 
at which stresses will be calculated and printed. 

NSTEPS - Number of load steps to be used in the solution 

process. A larger number of steps requires less 
iterations per step to reach convergence. The 
number of load steps which will result in minimum 
computational time is problem dependent. For 
efficiency, NSTEPS should be unity when only 
linear materials are used. 

MESH - Defines the refinement of a mesh of subtriangles 

used for numerical integration of stresses over a 
triangular region of a SIP or filler bar. Each 
side of the triangular region is divided into 
"MESH" equal intervals in forming the subtriangles. 
A linear distribution of stress is assumed over 
each subtriangle in the mesh. Therefore, greater 
accuracy but longer computation time results from 
a finer mesh. An example of the subtriangle 
pattern for MESH =4 is shown below. 
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TEST - The convergence test parameter for the iterative 

solution procedure. The iterations are terminated 
when the difference between each component 
(Fx» Fy, F2, Mx, My, M^) of the externally applied 
loads and pad material reactions is less than the 
value of TEST. 

NCON - Number of intervals along both axes of a rectangular 
grid over each pad at which stresses will be 
printed. Up to 34 intervals may be specified. 

If NCON is specified as zero, these stress dis- 
tributions will not be printed. The boundary of 
a rectangular grid is calculated from the minimum 
and maximum local x and y coordinate values of 
the corner points of the pad. If the pad itself 
is not rectangular, values on the grid outside 
the pad boundary will be printed, but should not 
be mistakenly interpreted as SIP stresses. Only 
values inside the pad boundary are SIP stresses. 

Command: "NONSTRESS" 

Purpose: The "NONSTRESS" command causes execution of the nonlinear 

analysis procedure. 

Format: "NONSTRESS" LI L2 ... / IMPl IMP2 ... / 

Description: The analysis will include the load case combination 

LI + L2 + ... and imperfections IMP1+IMP2+ ... This com- 
mand may be executed multiple times during a single run for 
different load conditions and imperfection combinations. 
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EXAMPLE ANALYSIS 


This section contains input data and corresponding calculated results for 
the analysis of an example tile and is designed to illustrate use of all input 
commands. This example tile is representative of an acreage tile having a 
6-inch square planform and a thickness of 2 inches. A graphical representation 
of the tile geometry is shown in figure 2. The tile could have been modeled 
using a single block element but is divided into a block, a wedge, and three 
tetrahedrons to illustrate use of various geometry input commands. The tile is 
mounted on a single SIP surface with corner nodes numbered 5, 6, 7, and 8. The 
SIP is surrounded by a filler bar with a width of 1/2 inch. 

The loads applied to the tile include a concentrated force and moment, a 
pressure gradient on the outer surface, and a statically equivalent acceleration 
or g-load. These loads are set up to be symmetric about the diagonal. Nodes 1 
and 3, of the tile. A nonlinear material curve is input for the SIP in the 
normal or through-the-thickness direction and linear properties are specified in 
shear. Filler bar properties are defined on a separate data file as illustrated 
in appendix A. A nonlinear analysis is performed for a combination of all four 
load cases and the sum of a cylindrical and double cosine imperfection. The 
input commands for this example tile analysis follows: 
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"TITLE" 

ACREAGE TILE EXAMPLE 

"CASEID" 1U222S33 1.05 635.3 3.63 

"COMMENT" 3 


THIS example 18 DESIGNED TO ILLUSTRATE 
USE OF ALL available INPUT COMMANDS. 
PREPARED BY GARY L. GILES 3«23*S1. 


"NODE" 

1 

0.0 

0.0 

0.0 



"NODE" 

2 

0.0 

6.0 

0.0 



"NODE" 

3 

6.0 

6.0 

0,0 



"NODE" 

H 

6.0 

0.0 

0.0 



"CORNER" 

ti 

1 

2 

5 0.5 

0.5 


"CORNER" 

1 

2 

3 

6 0.5 

0,5 


"CORNER" 

2 

3 

4 

7 0,5 

0,5 


"CORNER" 

3 

4 

1 

8 0,5 

0,5 


"NORMAL" 

u 

1 

2 

9 2.0 



" normal " 

1 

2 

3 

10 2,0 



"NORMAL" 

2 

3 

4 

11 2,0 



"NORMAL" 

3 

4 

1 

12 2,0 



"POINT" 

1 

4 

1 

13 3,0 



"POINT" 

9 

12 

9 

14 3.0 



"POINT" 

10 

11 

10 

15 3.0 



"PLAN£*LINE" 

13 

14 

15 

2 3 

16 


"LOCAL" 

1 

2 

3 




"BLOCN" 

1 

2 

16 

13 9 

10 15 14 

" WEDGE " 

3 

4 

U 

16 13 

15 / . 

005208 

"TETRA" 

11 

13 

14 

15 / .005206 


"TETRA" 

4 

11 

14 

12 / .005206 


"TETRA" 

4 

11 

14 

13 / .005206 


"FACE" 

9 

10 

11 

12 / .011 

.061 

17 

"FACE" 

1 

2 

10 

9 / ,011 

,061 

16 

"FACE" 

2 

3 

11 

10 / ,011 

,061 

19 

" FACE " 

3 

4 

12 

11 / ,011 

.061 

20 

"FACE" 

4 

1 

9 

12 / .011 

.061 

21 

"CG» 

22 






" VIEW " 

1 

2 

3 




"FORCE" 

1 

17 

5.0 

5,0 

0,0 


"MOMENT" 

2 

1 

• 15.0 

15,0 

0,0 


"PRESSURE" 

3 

9 

1.5 

12 

1,0 

10 

"PRESSURE" 

3 

10 

1.0 

12 

1.0 

11 

"GLOAO" 

4 

0. 

0 

0.0 50,0 




.005206 


0 

5 
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"MATEHIAL** I "TABIE* .160 «3 

••09000 •12.00000 

•*06600 •10.00000 

•.06200 •6.00000 

••07500 *6.00000 

-.06625 •0.68000 

-.06355 -0.26000 

•.06066 *3.80000 

-.05713 *3.39000 

-.05305 -2.95000 

•.00860 -2,53000 

•.00268 -2.06000 

-.03590 *1.62000 

-.02963 -1.29000 

-.02095 •1,07000 

•.01970 •, 65000 

-.01600 •, 73000 

•.01310 -.61000 

-.01060 -, 53000 

•,00736 -.02000 

-.00066 -,32000 

-.00303 -.26000 

-.00106 -.16000 

-.00056 -.11000 

0.00000 0.00000 

.00600 .25000 

.01100 .50000 

,01500 ,83000 

.01600 1,12000 

.02200 1, 69000 

.02700 2,99000 

.03200 3.63000 

.03600 5,10000 

.09000 6,80000 

,09600 9,10000 

.05000 10,90000 

,05500 11,70000 

,06000 12.60000 

.06300 13.90000 

.06500 13.90000 

.06600 13.30000 

.06800 12,75000 

.07000 11,30000 

,06000 0,00000 

••MATEHIAL" 2 •NOHMAL" "FBN" "N16C" 1000 0 

"MATEHIAL" 3 "LINEAR" 9,0 ,16 

"IMPERFECT* 1357 0,02 12,0 0,0 

"IMPERFECT" 2957 -0,0075 12,0 2,0 

,"PA0" 1 "F8" 2012 3 9158765/12/ 

"PAD" 2 "SIP" 135678/12/ 

"CONTRUL" 2 2 2 15 ,001 10 

"NONSTRESS" 1 2 39/1 2 / 
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The printed output listing is generated on the file, TAPE6. Selected 
portions of this output listing for the example tile will be presented and 
described. Portions of the output which are self-explanatory have been omitted 
for brevity. Also, portions which contain the same type of information as the 
included portions but are for different load steps and material pads have been 
omitted. 

As input commands defining geometry, loads, and materials are read, they 
are printed in the output along with some self-explanatory, calculated Informa- 
tion such as element weight and center-of-gravity location. This portion of the 
listing corresponds to all commands in the example input listing except for the 
"NONSTRESS" command and is not included herein. 

The "NONSTRESS" command initiates execution of the nonlinear stress 
analysis of the model defined by all preceding commands. Values for each of 
the defined load cases, referenced to the tile center of gravity, are printed 
as shown on the following page. These load vectors correspond to F^, Fy, F^, 

^ The global coordinate system. Next, the loads and imperfections 
which are to be combined for a particular analysis are printed with the corre- 
sponding total applied tile loads. 
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NQNSTRESS 

RESOLUTION OF EXTERNAL LOADS 

TU TILE C,G, 

a, 720 

-,000 

\ 

5,000 

5,000 

0,000 -«,720 

2 

0,000 

0,000 

0,000 -15,000 

15,000 

0,000 

3 

0,000 

0,000 

30,000 -9,000 

9,000 

0,000 


0,000 

0,000 

21,567 0,000 

0,000 

0,000 


ACREAGE TILE EXAMPLE 
nonlinear stress analysis • NASArLRC 

LOADS COMBINED I 3 3 A 

SELECTED IMPERFECTIONS# i, i$ 

EXTERNAL LOADS ■ .SOOOOE+Ol ,50000E+01 ,S7567E+0a ■,287a0Et02 ,2072OEfO2 •,7105«E»15 


CJ 

LD 



During the solution procedure, the total external tile load can be 
applied as a sequence of equal load increments or steps. These load steps are 
defined in the program by multiplying the total load, shown on the previous 
page, by a factor, lambda. An analysis is performed for the unloaded tile, 
lambda = 0.0, and in equal steps up to the total load, lambda =1.0. A con- 
verged solution is obtained at each load step. At each load step or value of 
lambda, the solution iteration history can be printed if selected by the value 
IPRT on the "CONTROL" command. The iteration history contains the values of 
displacements and rotations of the tile center of gravity (QCG) and correspond- 
ing values of the difference between the applied forces and pad material 
reactions (DELF). The iteration histories for the acreage tile example at the 
initial and final load steps are shown on the following page. The initial set 
of iterations, which are numbered one, for lambda equal zero correspond to 
reaching equilibrium with no applied loads or imperfections. This resulting 
equilibrium position is solely a function of the material stress displacement 
curves and allows curves which do not pass through the origin to be used. This 
position is taken to be the reference location from which all displacement 
values are calculated herein. The remaining iterations at lambda equal zero 
seek equilibrium for the specified substrate deflections. All remaining 
iteration histories are for applied loads. 
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ACREAGE TILE EXAMPLE 


LAMIjDAs 0,000 


ITER« 

1 

GCGp 

0, 

0, 

0, 

0 , 

0 , 

0, 



DELFo 

0, 

0, 

,15987E-13 

«,40390E«26 

,55538E-26 

0. 

ITER« 

1 

QCGb 

.12025E-29 

.U350E-29 

.30161E-17 

•,10748E-29 

,U766E-29 

•.31862E-57 



OELFa 

0, 

0, 

,15987E-13 

-,4O390E*26 

.55536E-28 

0. 

ITERii 

1 

QCGa 

,l2a25E*29 

.U3S0E-29 

.30161E.17 

•,10748E«29 

.U766E-29 

•,3l862E-57 



UElFa 

•,3l529E-«0 

-,31529E«40 

,29339E+02 

,86509E400 

•,66509E400 

•46940E-S4 

ITERt 

2 

QCGa 

-,85459E«00 

•,eS459£*04 

,19333£«01 

,80927E-04 

.,eo92?E-04 

•.246S2E-31 



CELFa 

,ta296E*l4 

,30493E*i4 

.U918E + 01 

-,32297E+00 

.32297E+00 

•,27768E-27 

ITER« 

3 

QCGa 

•13282E*05 

,13282E«05 

,20186E-01 

-.12578E.05 

,12578E*05 

•,17256E«30 



DELFf 

-,77186E-15 

-,71468E-15 

.71863E-01 

•.U282E+00 

,U282Ef00 

•i08S5E«26 

ITERa 

<i 

QCGa 

,4086eE»Q4 

,40e6BE*04 

,20248E-01 

•,38700E»04 

,38700E-04 

•12326E-31 



CELFa 

-,30a93E-l5 

,30493E«15 

-,3276lE-02 

,30020E»02 

•,30020E«02 

•,10729E»27 

ITERS 

5 

QCGa 

,39820E«oa 

,39820E*04 

,20246E*01 

•,37708E*Oa 

.37708E-04 

••44296E-31 

• 

0 


CELFa 

•,30a9JE-lS 

,7892lE»4l 

-,27759E-05 

.,35099E-05 

,35099E-05 

r56798E-26 


LAMUDA* 1,000 


ITERa 

1 

QCGa 

,54960E«02 

.54960E«i02 

,3480«E-0l -,35210E*02 

,35210E-02 

•17661E*16 



CELFa 

,25000E401 

,25000E*01 

• 26784Et02 -.HSEOEfOE 

•14360E402 

•56843E013 

ITERS 

2 

QCGa 

,«6570E*02 

,86S70E»02 

,43678E-0l •,48309E»02 

,48309E«02 

.45023E-16 



CELFa 

-,25580£«12 

•,22737E-12 

•.99901E+01 .62582E+01 

-•62S82E401 

,24158E-12 

ITERa 

3 

QCGa 

,83247E-02 

,832a7£*02 

,4l957E-0l -,45162E*02 

•45162E-02 

.86253E-16 



CELFa 

-,17053E-12 

-,170S3E.12 

•,41643E+00 ,20044E+00 

•,2O044Ef00 

•,75318E-12 

iter* 

4 

QCGa 

.83165E-02 

,83l65E-02 

,41873E-01 -,4S084E»02 

,45084E*02 

•.42289E-16 



CELFa 

,5b8a3E-13 

-.28422E-13 

-,87313E-03 .20720E-03 

•,20720E*03 

,78160E»12 



Stress information is always printed for the total applied load, lambda 
equal 1.0, and can be optionally selected at each of the load steps. This 
stress output is printed separately for each pad and can be identified by the 
user defined pad number. The output for the SIP, pad 2, for the example tile 
is given on the following page. The material curve indicators used for normal 
and shear calculations are printed with imperfection indicators and the applied 
load cases. Forces in the pad local coordinate system which are equivalent to 
the pad stresses are then output. 

Next, a table of stress related data which are computed during the analysis 
is printed. The information is printed at the corners of the pad, identified 
by node number, and any intermediate locations on the pad edge which were 
selected by NDIV on the "CONTROL" command. 

The coordinates of the locations where information is to be output are 
given by the values XX and YY in the pad local coordinate system. Displace- 
ments at the inner surface of the tile (XDISP, YDISP, and ZDISP) are given in 
the pad local coordinate system relative to the equilibrium condition with no 
imperfections and no externally applied tile loads. Shear stresses, TAUXZ and 
TAUYZ, are given in the local X-Z and Y-Z planes, respectively. The normal or 
through-the-thickness stress in the pad is denoted as SIGZ. The total substrate 
deflection, WIMP, is the sum of the "IMPERFECT" contributions and the total 
normal displacement of the pad material relative to the origin of the stress- 
displacement curve is given by WEXT. The values DTAUXZ, DTAUYZ, and DSIGZ are 
the slopes or derivatives of the corresponding stress-displacement curves for 
the displacement value at that particular pad location. 
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ACREAGE TILE EXAMPLE 


PAD NOMBtfi 2 

material * NORMA 
material - SHcAR 
IM.pEkMECTICNS « 
loads COMblNECs 

L s 1 

* 3 

1# 

»» 

*SIF * 

St 

St It Ht 




RESOLOIICN OF FORCES ON 
LOCAL COCRCINATES 

PAD* .50000E + 01 .SOOOOE-FOl .57638E + 02 ,1l025E-» 

03 110255403 ,176335 

node 

XX 

TV 

XOISP YDISP ZDI3f> 

RIMP WEXT 

TAUXZ 

DTAUXZ 

TAUYZ 

OTAUYZ 

SIGZ 

OSIGZ 

s 

0.000. 

0,000 

.35556E-02 .35556E-02 .60015E-01 

.12500E-01 ,51915E-0i 

.20000E400 

.56250E402 

.20000E400 

.562505402 

.10698E402 

.260005403 


0.000 

2,500 

.3S556E-02 ,35SS6E"02 ,S3ia4E«01 

,2R166E»01 ,269765-01 

,20000E«00 

,562505402 

,200005400 

.56250E402 

.29109E401 

.23600E403 

b 

0.-000 

5,000 

,3S556E-02 .3S556E-02 .01673E-01 

,666506-02 ,33206E-01 

,200005400 

.562505402 

,200005400 

,562505402 

,40739E401 

,367505403 


2.500 

5.000 

,3S556E-02 .355S6E-02 .30602E-01 

,2097AE-01 .96279E-02 

,200005400 

.56250E402 

,200005400 

.56250E402 

.O3139E4O0 

,500005402 

7 

■5.000 

S.ooo 

•35556E-02 .35556E-02 .19331E-01 i 

.19160E-01 .166O5E-03 

.200005400 

,562505402 

.20000E400 

.562505402 

,691885-02 

,016675402 


•s.ooo 

2,500 

•35556E-02 .355S6E-02 .30602E-01 

.20970E-01 ,962796-02 

.200005400 

.S6250E402 

,200005400 

.562505402 

,431395400 

,500005402 

a 

•5.000 

0,000 

•35SS6E-02 ,35S56E-02 ,01673E-01 

.20000E400 

•20000E400 

,407395401 




.66650E-02 ,33208E-01 

.562505402 

.S6250E402 

.367505403 


2,500 

0,000 

•35556E-02 ,35SS6E-02 .S3144E-01 

, .24166E-01 ,2697eE-01 

,200005400 

,562505402 

.200005400 

.562505402 

,291095401 

,238005403 

5 

0.000 

0,000 

• 35556E-02 ,35S56E-02 ,6<t415E-01 

•12500E-01 ,S1915E-01 

,200005400 

,562505402 

.20000E400 

.562505402 

,106985402 

,260005403 


0,000 

0,000 

• 3S5S6E-02 ,3S556E-02 ,64<I15E-01 

• 12SO0E-01 .51915E-01' 

.20000E400 

•S6250E402 

.20000E400 

.S6250E402 

,106985402 

,260005401 



Displacement information is printed next for the total loading condition 
as shown on the following page for the acreage tile example. Displacement 
vectors, three translations, and three rotations are given at the tile center 
of gravity for the reference equilibrium position (lambda = 0 and no imper- 
fections) and for the total applied load condition (lambda = 1.0 and including 
imperfections). A table of displacements, DELX, DELY, and DELZ, are then 
printed for each node point. Also included in the table are the user defined 
coordinates of the nodes. The displacement values correspond to the difference 
between the reference and final equilibrium positions. All values are given 
in the global coordinate system. 


44 



acreage tile exa^<ple 


FINAL NODAL DISPLACEMENTS IN 
QCGREF AT LAMBDAOOfANO NO IHpERFECTIONSs 


. 12625 E 

-29 

11350 E -29 

. 30161 E -17 

10768 E -29 

AT LAMBOAsl , including ANY 

imperfections # 

.831656 

>02 • 

83165 E -02 

. 61873 E -01 

•. 65085 E -02 

NODE 

X 

Y 

z 

DELX 

1 

0.0000 

0.0000 

0.0000 

.0036 

2 

0.0000 

6.0000 

0,0000 

.0036 

3 

6.0000 

6,0000 

0,0000 

.0036 

6 

6.0000 

0,0000 

0,0000 

.0036 

5 

,5000 

.5000 

0,0000 

,0036 

6 

.5000 

5,5000 

0,0000 

,0036 

7 

5,5000 

5,5000 

0,0000 

.0036 

8 

5,5000 

,5000 

0,0000 

.0036 

9 

0.0000 

0,0000 

2,0000 

,0126 

10 

0.0000 

6,0000 

2,0000 

.0126 

1 1 

6.0000 

6,0000 

2,0000 

.0126 

12 

6.0000 

0,0000 

2,0000 

.0126 

13 

3,0000 

0,0000 

0,0000 

.0036 

16 

3.0000 

0,0000 

2,0000 

.0126 

15 

3.0000 

6,0000 

2,0000 

.0126 

16 

3.0000 

6,0000 

0,0000 

.0036 

17 

3.0000 

3,0000 

2.0000 

.0126 

18 

0.0000 

3,0000 

1,0000 

.0081 

19 

3.0000 

6,0000 

1,0000 

.0081 

20 

6.0000 

3,0000 

1,0000 

.0081 

21 

3,0000 

0.0000 

1,0000 

,0081 

22 

3,0000 

3.0000 

1,0560 

,0083 


on 


SYSTEM 


, U 766 E . a «> -. 31862 E -57 
, iJ 5085 E -02 . 9 U 0 aE - l 6 


OELY 

.0036 

• 0036 
.0036 
.0 036 
.0036 

• 0036 
.0036 
.0036 
.0126 
.0126 
.0126 
.0126 
.0036 
.0126 
.0126 
.0036 
.0126 
.0061 
.0061 
.0061 
.0081 
.0083 


DELZ 

.0669 

.0619 

.0168 

.0619 

.0666 

.0619 

.0193 

.0619 

.0689 

.0619 

.0168 

.0619 

.0556 

.0556 

.0283 

.0283 

.0619 

.0556 

.0283 

.0283 

.0556 

.0619 



A nonzero value of NCON on the "CONTROL" command is used to request print- 
ing of imperfection and stress values at points on a grid over each SIP pad. 

The boundary of this rectangular grid is calculated from the minimum and maximum 
X and Y values of corner nodes in the pad local coordinate system. If the pad 
is not rectangular, values on the grid outside the boundary will be printed but 
should not be mistakenly interpreted as SIP stresses. 

The output grids of imperfection and stress data for the acreage tile 
example are shown on the following page. The origin of the local pad coordinate 
system corresponds to the value in the lower, left corner of the table with the 
local y-axis positioned vertically and the local x-axis positioned horizontally. 
The maximum and minimum stress values are calculated and printed on separate 
lines. For convenience of use at an interactive terminal, the stress grid and 
associated maximum and minimum stresses for each pad are output on a file 
named TAPEl . 

When multiple tiles and/or multiple load conditions are analyzed during a 
single run, the information of primary interest is the maximum stress values 
that are calculated. Such information is output on a file named TAPE2 which 
contains a table of values with each line having information from the CASEID 
command, a load condition identifier, and a maximum stress value. Each line 
of the table is written as a tile part number, the load identifier, Mach number, 
dynamic pressure, angle of attack, and maximum stress as shown below for the 
acreage tile example. 

111222333 I i,05 635,3 3,83 10,8979 
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ACREAGE tile example 


imperfection grid sip 2 








LOADS 

COMBINED* 

1» 

20 

009 

,008 

,010 

,012 

,016 

.021 

SELECTED 
.025 .026 

imperfections* 

.026 .023 

1 0 
,019 

oia 

,011 

,010 

.011 

,019 

,016 

.023 

,026 

,027 

,026 

.023 

016 

.019 

,012 

• Oil 

,012 

.015 

,019 

,023 

,026 

.027 

,026 

oao 

,019 

,016 

• 013 

.012 

.013 

,015 

.019 

,023 

• 026 

,026 

023 

,023 

,021 

,017 

,019 

,013 

.013 

,015 

• 019 

.023 

.025 

029 

.026 

,025 

,022 

.016 

,019 

,013 

,013 

,015 

,016 

.021 

023 

,026 

,027 

.025 

,022 

,016 

,019 

.012 

,012 

,019 

.016 

020 

.029 

,027 

,027 

.025 

.022 

,017 

,013 

,011 

• Oil 

.012 

016 

,021 

,025 

,027 

,027 

.025 

.021 

,016 

.012 

,010 

,010 

019 

.017 

,021 

,029 

.026 

.026 

• 023 

,019 

.019 

• Oil 

,006 

013 

,019 

,016 

,020 

.023 

.029 

,023 

,020 

,016 

,012 

,009 


STRESS GRID SIP 2 


9,07 

3.92 

2.62 

1.^5 

.96 

.93 

3,67 

3,95 

3,05 

2.33 

1.53 

,81 

3,19 

3,12 

3,07 

2,76 

2.15 

1,90 

2.75 

2,60 

2.71 

2.76 

2,56 

2.05 

2,59 

2.25 

2,26 

2.93 

2.61 

2,95 

2,91 

2,16 

1,98 

2,06 

2.33 

2.55 

3.62 

2,96 

1,96 

1.69 

2.02 

2.33 

5.89 

3,96 

2.32 

1,88 

1,69 

2,06 

6,15 

5,59 

3,36 

2.32 

1,96 

1,96 

9,99 

6,09 

5,59 

3,96 

2.96 

2,16 

10,90 

9,99 

8,15 

5.69 

3.62 

2,91 


MAXIMUM STRESS ■ 10,90 
MINIMUM STRESS • -.aS 


LOADS COMBINEDr 1, 2, J, 

SELECTED IMPERFECTIONS* i, 2 , 


.15 

• ,07 

•,21 

-,18 

.01 

.35 

.10 

-.17 

-.23 

-.18 

.73 

.31 

,08 

•.17 

-.21 

1,33 

.70 

.31 

• 10 

-.07 

2.01 

1,33 

.73 

.35 

• 15 

2.95 

2,05 

1.90 

.61 

.93 

2.61 

2,56 

2,15 

1.53 

,96 

2.93 

2,78 

2,76 

2.33 

1.75 

2,26 

2.71 

3,07 

3.05 

2.62 

2,25 

2,60 

3.12 

3,95 

3,92 

2,59 

2.75 

3,19 

3,67 

9,07 




END OF RUN 



The tile part number, Mach number, dynamic pressure, and angle of attack are 
values input on the CASEID command. The load identifier is the first load case 
given on the NONSTRESS command and the last value is the corresponding calcu- 
lated maximum stress. 

The printed output described in this section is repeated for each NONSTRESS 
command which is input. Termination of the run is indicated by the END OF RUN 
message as shown on the bottom of the preceding page. 

CONCLUDING REMARKS 

User documentation is presented for a computer program which considers the 
nonlinear properties of the strain isolator pad (SIP) in the static stress 
analysis of the shuttle thermal protection system (TPS). This program is 
generalized to handle an arbitrary SIP footprint including cutouts for instru- 
mentation and filler bar. Multiple SIP surfaces can be defined to model tiles 
in unique locations such as leading edges, intersections, and penetrations. 

The nonlinearity of the SIP is characterized by experimental stress-displacement 
data for both normal and shear behavior. Stresses in the SIP are calculated 
using a Newton iteration procedure to determine the six rigid body displacements 
of the tile which develop reaction forces in the SIP to equilibrate the 
externally applied loads. 

This user documentation gives an overview of the analysis capabilities, a 
detailed description of required input data, and an example to illustrate use 
of the program. 


48 



APPENDIX A 


MATERIAL CURVE DATA FILE 

The pad material curves used in the analysis are a function of material 
thickness, orientation (normal or shear directions), and the history of loading 
the material has experienced. Thus, a large number of material curves may be 
required for use in an analysis effort. The SPAR data handling utilities, 
reference 4, can be used by the analysis program for convenient bookkeeping, 
storage, and retrieval of these material curves. The curves are stored on a 
file called SPARLA and each curve is referred to by a unique, four word name 
assigned by the user. An example of input data required to establish or add 
curves to the SPARLA data file is given in this appendix. 

The input statements are free-field format and a detailed description of 
the statements used is given in reference 7. The example input which follows is 
for the filler bar curve which was used for the acreage tile example presented 
in the main body of this report. 
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[XQT AUS 
TABLE(NI=2,NJ=29) 
FBN N16C 1000 0 
1 = 1,2 


J=l,29 

0.16 

0.0 

-.07876 

-9.95000 

-.07640 

-8.96000 

-.07490 

-8.36000 

-.07299 

-7.69000 

-.07106 

-7.10000 

-.06870 

-6.46000 

-.06677 

-5.95000 

-.06386 

-5.32000 

-.06120 

-4.87000 

-.05847 

-4.35000 

-.05537 

-3.88000 

-.05242 

-3.49000 

-.04909 

-3.09000 

-.04535 

-2.72000 

-.04106 

-2.31000 

-.03732 

-2.01000 

-.03284 

-1.70000 

-.02850 

-1 .41000 

-.02423 

-1.16000 

-.02053 

-.95000 

-.01722 

-.79000 

-.01343 

-.61000 

-.01112 

-.50000 

-.00791 

-.37000 

-.00453 

-.21000 

-.00177 

-.09000 

0.00000 

0.00000 

.01600 

0.00000 

TABLE{ 

) 


[XQT 

DCU 

TOC 1 


PRINT 

1 

[XQT 

EXIT 


Each curve is given a four word name; the first two are alphanumeric and 
the last two are integers. The name of the curve in this example is 
FBN N16C 1000 0 as shown in the third statement. The values NI and I 
refer to the number of columns in the table which is two. Values NJ and J 
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refer to the number of rows in the table which is twenty-nine for the example 
shown. The first row of the table contains the thickness of the pad and a 
dummy value (0.16 and 0.0 in the example). The remaining twenty-eight rows 
are points on the stress-displacement curve with the displacement value given 
first followed by the stress value. Such tables of data can be repeated for 
input of multiple curves as indicated. The last four statements shown in the 
example request the information on the data file to be printed for verification 
and termination of the run. 
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(a) TPS arrangement. 



(b) Tile modeling. 


(c) Nonlinear modeling. 
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Figure 1.- Analytical model of thermal protection system. 
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Figure 2.- Geometric model of acreage tile example. 
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